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Perfect and broadband acoustic absorption by critical coupling 
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We experimentally and analytically report broadband and narrowband perfect absorption in two 
different acoustic waveguide-resonator geometries by the mechanism of critical coupling. In the 
first geometry the resonator (a Helmholtz resonator) is side-loaded to the waveguide and it has 
a moderate quality factor. In the second geometry the resonator (a viscoelastic porous plate) is 
in-line loaded and it contains two resonant modes with low quality factor. The interplay between 
the energy leakage of the resonant modes into the waveguide and the inherent losses of the system 
reveals a perfect and a broadband nearly perfect absorption. The results shown in this work can 
motivate relevant research for the design of broadband perfect absorbers in other domains of wave 
physics. 
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Structures composed of waveguides and resonators 
have been exploited for fundamental studies as well as for 
several applications in different branches of science and 
technology [1] . In these systems the balance between the 
leakage rate of energy out of the resonator (resonator- 
waveguide coupling) and the inherent losses of the res¬ 
onators has been shown of fundamental relevance for 
their transmission and reflection properties [m. When 
the leakage and the inherent losses are well balanced, the 
critical coupling condition is fulfilled, and then, a perfect 
destructive interference between the transmitted and the 
internal fields leads to maximum absorption at the res¬ 
onance frequency BIS]. Therefore the critical coupling 
paved the way to produce perfect absorption (PA) of the 
incoming waves. In particular, the coherent perfect ab¬ 
sorption [6] , being the time-reversed counterpart to laser 
emission [7], makes use of the critical coupling mecha¬ 
nism to perfectly absorb the two port coherent excita¬ 
tion in a lossy medium. Recently, it has been shown that 
PA can be obtained by a single port excitation in a two 
port structure [8]. This is accomplished by the degen¬ 
erate critical coupling condition, i.e., if two degenerated 
resonances are critically coupled and have opposite sym¬ 
metry. 

PA is of particular interest for many applications such 
as energy conversion 0 uni, time-reversal technology 
[TTJ [12] , coherent perfect absorbers laizi or soundproof¬ 
ing [131 IH] among others. In optics, several configu¬ 
rations using Bragg reflectors [Tbj, Fabry-Perot cavities 
with metamaterial mirrors m, layered media with Kerr 
nonlinearity m or graphene-based hyperbolic metama¬ 
terials [IHI have been recently used to obtain nearly PA 
of the incident radiation by critical coupling. Directional 
perfect absorbers have been also designed using deep sub¬ 
wavelength low-permittivity films m- Moreover PA has 
been obtained in plasmonic planar structures via the crit¬ 
ical coupling of the surface plasmon modes [20] , while this 
has been also used for heat generation in plasmonic meta¬ 
materials [21]. In acoustics, coherent perfect absorption 


has been studied numerically by controlling the coher¬ 
ence of the input waves in a two port system showing 
a dynamically tuned absorption coefficient on the target 
material up to unity [22]. 

In this study, we exploit the critical coupling to exper¬ 
imentally and theoretically show narrowband and broad¬ 
band PA in two simple configurations: a side or an in-line 
resonant element coupled with a waveguide backed by a 
rigid end. The frequency range considered here is well 
below the first cut-off frequency of the higher propaga¬ 
tive modes in the waveguide, therefore the problem is 
considered one-dimensional. Both configurations can be 
considered as an asymmetric Fabry-Perot cavity of length 
L with two different mirrors, i.e., the resonant element 
and the rigid backing [as schematically shown in Fig. 
I^a)]. The absorption of this system can be expressed 
as a = I — |rp, where r is the complex reflection co¬ 
efficient obtained from the standard three-medium layer 
Fresnel equation [23] . 

I — 

where the time harmonic convention is vr and 

Vt are the reflection coefficients of the resonant element 
and of the termination, respectively (in our case, = I), 
and k = uj/c is the wave number in the cavity of length 
L. Gonsidering losses in this configuration, the PA is 
fulfilled when the reflection coefficient is zero, i.e., when 
the superposition of the multiple reflections in the cavity 
[second term in Eq. Q] destructively interferes with the 
direct reflection from the resonant element (first term in 
Eq. 0). 

We use Helmholtz resonators (HRs) to build the side 
resonant coupling geometry [see Fig. [^b)], and viscoelas¬ 
tic porous plates to build the in-line resonant coupling 
geometry [see Fig.j^c)]. We choose these two systems to 
analyze two different situations: a weakly lossy single res¬ 
onant element with moderate quality factor {Q~^ ^ O.I, 
in the case of HRs) and a largely lossy double resonant 
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FIG. 1. (Color online) (a) Asymmetric Fabry-Perot resonator 
made of a resonant element (red point) and a rigid backing 
at distance L from the resonator, (b) and (c) respectively 
show the side and in-line resonant coupling geometries backed 
by a cavity of length L. (b) experimental set-up using an 
impedance sensor (IS) as described in Ref. [24]. The HR is 
composed of a neck of length Ln — 2 cm with radius Rn — 1 
cm, a cavity with tunable length, Lr^ and radius Rr = 2.15 
cm. The waveguide has a radius R = 2.5 cm. (c) experimental 
set-up based on the ISO-10534-2 m- The viscoelastic porous 
plate of radius Rp = 2.2 cm is clamped in a tube of R = 2.2 
cm. 


element with low quality factor {Q~^ — 1, in the case 
of the viscoelastic porous plate). We show that PA can 
be achieved though the critical coupling of the resonant 
modes of the system. In the largely lossy case, we take 
advantage of the losses to overlap PA peaks with other 
nearly PA peaks and thus to obtain a broadband nearly 
PA. The two configurations have been experimentally 
analyzed using standard impedance tubes showing very 
good agreement with the theoretical predictions (see cap¬ 
tion of Fig. for references and details of the set-ups). 

To understand the origin of the PA in our systems, 
we analyze the reflection coefficient r evaluated in the 
complex frequency plane, i.e., substituting uj = ujr + lUJi 
in k. This kind of map has been revealed extremely in¬ 
formative for understanding the PA process [T9|. With¬ 
out losses, one finds pairs of poles (at complex frequency 
ujp) and zeros (at complex frequency uJz = ojp) of r 
symmetrically distributed around the real frequency axis 
m- The zeros (poles) are in the positive (negative) half 


imaginary frequency plane EH- When the losses are in¬ 
troduced into the system, the zeros and the poles are in 
general down-shifted. In the case where the zero is on 
the real frequency axis {uJz = xr-\-iO = 27r/), r(ujz) = 0 
so the PA condition is fulfilled and the PA is obtained. 

We start by the analysis of the side resonant coupling 
geometry made of a HR with a variable cavity length, 
Lr^ attached to a fixed backing cavity of L = 15 cm [see 
Fig. Bb)]. This system presents viscothermal losses at 
the walls of the waveguide and of the resonator [24l |28] . 
By changing Lr from 0 to 15 cm, i.e., by changing the 
resonant frequency of the HR, we study the trajectory 
of the zero of the reflection coefficient in the complex 
frequency plane for both the lossless (black dashed line) 
and the lossy (black continuous line) cases in Fig. [^a). 
On one hand, as shown by the arrows over the trajec¬ 
tory, the zero moves to lower real frequencies as Lr in¬ 
creases. On the other hand, the presence of the losses 
down-shifts the trajectory of the zero with respect to the 
lossless case and thus, the PA condition can be satisfied. 
In this case, there are two crossing points with the real 
frequency axis, i.e., two different configurations accom¬ 
plishing the PA condition. These points correspond to 
(Lr,/) = (8.3 cm, 484.5 Hz) and (Lr,/) = (3.9 cm, 647 
Hz). The blue and red areas in Figure |^a) respectively 
show the zero and the pole of the reflection coefficient 
[Eq. Q] in the complex frequency plane for the config¬ 
uration (Lr,/) = (8.3 cm,484.5 Hz). This case corre¬ 
sponds to a configuration at which the zero is on the 
real frequency axis, therefore the PA is physically ob¬ 
servable. Similarly, a complex map for the configuration 
(Lr^/) = (3.9 cm, 647 Hz) with the zero in the real fre¬ 
quency axis can also be obtained E3- 

We now theoretically and experimentally evaluate the 
absorption coefficient a for these two configurations cor¬ 
responding to the crossing points with the real frequency 
axis. We find analytically PA (a = 1) for the two above 
mentioned configurations at 484.5 Hz and 647 Hz as 
shown in Fig. ib), in agreement with the crossing of 
the trajectories with the real frequency axis represented 
in Fig. I^a). Experiments show also very good agree¬ 
ment with the theoretical predictions, producing 100% 
of absorption for these configurations at the correspond¬ 
ing frequencies of the PA condition. 

Einally, the theoretical and experimental dependencies 
of a on both Lr and / are plotted in Eigs. [^c) and|^d) 
respectively. Several configurations, around the PA con¬ 
dition, are able to produce high absorption peaks in this 
system. In these configurations, the zero of the reflec¬ 
tion coefficient is located close to the real frequency axis 
and influences the absorption coefficient that is physically 
observable. To characterize this high absorption, we de¬ 
fine the region of nearly PA using the isoline a = 0.9 
(|rp = 0.1) around the zero in the complex frequency 
plane. If this isoline crosses the real frequency axis, a 
nearly PA peak with absorption a > 0.9 becomes ob- 
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FIG. 2. (Color online) Side resonant geometry, (a) represents the complex frequency map of log|r| for the side resonant 
configuration with Lr = 8.3 cm. Black dashed (Continuous) line represents the trajectory of the zero of \r\ for the lossless 
(lossy) case in the complex plane as Lr is increased (sense of the increasing shown by arrows over the lines). White continuous 
line shows the isoline a = 0.9. Inset shows a zoom of the complex frequency map. (b) Green dashed (green open squares) and 
red continuous lines (red open circles) show the theoretical (experimental) absorption coefficient, a = 1 — |r|^, for Lr = 3.9 
cm and Lr = 8.3 cm respectively (horizontal lines in (c) and (d)). (c) and (d) show the theoretical and experimental maps 
respectively for the dependence of a on both the Lr and the frequency, /. 


servable m- This is translated in Figs. ic) and [^d) 
by the two regions with high absorption around the PA 
condition. 

The study now will be focused on the in-line resonant 
coupling geometry [Fig. [^c)]. This system presents 
larger inherent losses than the previous case, and also, 
it has two resonances in the studied frequency range. 
For that, a viscoelastic porous plate with losses com¬ 
ing from viscoelasticity and porosity is used as shown 
in Fig. [^c). The elastic properties are modelled by a 
Voigt rheological model, in which the Young’s modulus 
is complex with a frequency dependent imaginary part, 
E = Eq ^ irjuj {Eq = 220 kPa and = 7 Pa s). Due to 
the porosity, the mass density is also complex with the 
form p = pp — %xl{pp = 28 kg/m^ and x = 350 kg 
m3 gi/2^ both cases, p and x are extracted from the 
experimental results m- Recently, a similar structure 
has been analyzed m showing that a small amount of 
losses can lead to the design of subwavelength impedance 
matched systems. In our work, the inherent losses are 
substantially larger and the system contains two resonant 
modes. Therefore, the effects of each mode can overlap 
to produce a broadband nearly PA. 

In Fig. [^a), we analyze the trajectory of the zeros 
of Eq. 0 for both cases, the lossless and the lossy, as 
we change the cavity length L from 0 to 13 cm. As 
shown by the arrows over the trajectories, the two ze¬ 
ros move to lower real frequencies as the cavity length 
L increases. When the inherent losses are taken into ac¬ 
count, the trajectory of the zero of the first mode crosses 
the real frequency axis corresponding to the configura¬ 
tion (T,/) = (2.71 cm, 820 Hz). Moreover, the trajec¬ 
tory of the zero of the second resonance crosses the real 


frequency axis in two points, corresponding to the con¬ 
figurations (T,/) = (1 cm, 1703 Hz) and (T,/) = (1.37 
cm, 1508 Hz). The color map in Fig. [^a) depicts the 
evaluation in the complex frequency plane of Eq. 0 for 
the configuration (T,/) = (2.71 cm, 820 Hz). The zero 
of the reflection coefficient of the first resonant mode is 
on the real frequency axis, while for the second resonant 
mode it is far away from the real frequency axis. How¬ 
ever, the isoline a = 0.9 around the zero of the second 
resonant mode crosses the real frequency axis, denoting 
a substantial increase of the absorption in this frequency 
range. Eigurej^b) represents the absorption coefficient 
evaluated for the configurations having the zero of the 
reflection coefficient in the real frequency axis, i.e., with 
the PA condition activated. Eor the configuration with 
L = 2.71 cm (black line), the absorption peak due to the 
first mode is a signature of the PA (experimentally 100% 
of absorption), while the second one produces a nearly PA 
peak explained by the crossing of the isoline a = 0.9 with 
the real frequency axis [see Eig. |^a)]. Therefore the zero 
produced by the second mode, although being far away 
from the real frequency axis, influences the absorption co¬ 
efficient and overlaps with the PA peak of the first mode 
producing a broadband absorption peak with nearly PA. 
The other configurations, L = 1 cm and L = 1.37 cm, 
clearly show the PA peak produced by the critical cou¬ 
pling of the second resonant mode [see arrows in Eig. 

ia)]. 

Einally, in Eigs. [^c) and[^d), we respectively plot the 
theoretical and experimental dependency of a on both L 
and /. Two zones, corresponding to the first and sec¬ 
ond resonant modes, with very high absorption around 
the PA configurations are visible. PA peaks appear only 
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FIG. 3. (Color online) In-line resonant geometry, (a) represents the complex frequency map of log \r\ for the in-line resonant 
conhguration with L = 2.71 cm. Black dashed (Continuous) and red dashed (continuous) lines represent the trajectory of the 
zero of |r| for the first and second resonant modes for the lossless (lossy) case as Lr is increased (sense of the increasing shown 
by arrows over the lines). White continuous line shows the isoline a = 0.9. Inset shows a zoom of the complex frequency map. 
In (b) green dot-dashed (green open triangles), blue dashed (open blue squares) and black continuous (black open circles) lines 
show the analytical (experimental) absorption coefficients for the cases L = 1 cm, L = 1.37 cm and L = 2.71 cm (horizontal 
lines in (c) and (d) respectively), (c) and (d) show the theoretical and experimental maps for the dependence of a on both the 
length of the cavity and the frequency. 


for three particular configurations [horizontal lines, cor¬ 
responding to the results shown in Fig. ib)] , but the 
areas around them still present very high absorption. 
It is worth noting here that, in contrast to the case of 
the HRs with moderate quality factor, the viscoelastic 
porous plates have a low quality factor. Therefore, even 
for the configurations with the zero of the reflection coef¬ 
ficient out of the real frequency axis, the isoline a = 0.9 is 
wide enough to cross the real frequency axis and generate 
broad areas with nearly PA. 

In conclusion, we have exploited the critical coupling 
to experimentally and analytically show narrowband 
and broadband PA of acoustic waves in two waveguide- 
resonator geometries: the side and the in-line resonant 
coupling geometries. The complex frequency maps of 
the reflection coefficient and the trajectories of its zeros 
in the complex frequency plane have been revealed very 
insightful. The presence of losses produces a down-shift 
of the lossless trajectory of the zero inducing its crossing 
with the real frequency axis. The PA occurs when the 
zero of the reflection coefficient is on the real frequency 
axis. Nearly PA becomes observable when the isoline 
Q/ = 0.9 crosses the real frequency axis. We have shown 
here that, making use of the losses of the system, one can 
overlap PA and nearly PA peaks to produce broadband 
PA absorption. This could allow challenging applications 
and advances in different domains of wave physics. 
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